Electron gas at the interface between two antiferromagnetic insulating manganites 
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We study theoretically the magnetic and electric properties of the interface between two antiferro- 
magnetic and insulating manganites: Lao.sCao.sMnOa, a strong correlated insulator, and CaMnOs, 
a band-insulator. We find that a ferromagnetic and metallic electron gas is formed at the interface 
between the two layers. We confirm the metallic character of the interface by calculating the in- 
plane conductance. The possibility of increasing the electron gas density by selective doping is also 
discussed. 

PACS numbers: 75.47.Gk, 75.10.-b, 75.30.Kz, 75.50.Ee 



I. INTRODUCTION 

In recent years, atomic control has been gained in the 
growth of multilayers of thin oxide films with different 
electronic and magnetic properties.— Interestingly, 'elec- 
tronic reconstruction' (or redistribution of charges) at 
the interfaces of strongly correlated systems^ has led to 
the appearance of phases showing a very different be- 
havior from that in the bulk of the constituent materi- 
alsf^i^i^i^i^i^i^i^i^. For example, superlattices of thin lay- 
ers of LaTiOa and SrTiOa show metallic conductivity, al- 
though both materials are insulators in bulkii SrTiOa is a 
band insulator while LaTiOs is a Mott insulator and has 
one extra electron on the Ti. This extra charge spreads 
across the interface producing the metallic behavior.— 
Another interesting example is the formation of a high 
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FIG. 1: (Color online) 2-dimensional (xz plane) projection 
of the cubic lattice considered. The circles represent the Mn 
sites while the squares are the A-sites (Ca or La) shifted by 
(1,1,1)1 with respect to the Mn. a is the lattice parameter. 
Full squares represent the Lag's CaQi^O^" plane, with a charge 
density of -1-0.5 per A atom. Empty squares correspond to the 
Ca^'''0^~ plane, which is neutral. Between, the gray (light 
blue) squares correspond to the average -1-0.25. We consider 
periodic boundary conditions in all three directions. 



mobility electron gas at the interface between the two 
band insulators LaAlOa and SrTiOs.— This electron gas 
has been shown to have magnetic properties^ and to be- 
have as a two dimensional superconductor at low tem- 
peratures The electric field effect, which is the basis for 
semiconductor transistors, can also be used to tune the 
charge density at interfaces between strongly correlated 
oxides, changing their properties. 

The redistribution of charge in Mott-band insulator 
heterostructures and the formation of a metallic inter- 
face have been studied theoretically using different tech- 
niques including Hartree-Fock theory to account for the 
on-site correlations,^'^*' dynamical mean field theory,^^ 
and Lanczos diagonalization of a quasi-one-dimensional 
lattice.— The common observation is the formation of a 
metallic interface region up to three-unit cells wide. In- 
terfaces between two different Mott insulators have also 
been theoretically studied!^ a two-dimensional metal is 
formed with a density of states that can be controlled by 
a remote (a few unit-cells away) doping layer. 

Perovskite manganites of formula A^+A'^^^MnOa 
(with n the density of electrons) are strongly correlated 
magnetic materials that present a great variety of phases 
depending very strongly on the close competition be- 
tween different interactions: insulating behavior is fa- 
vored by electron-lattice interaction, Coulomb interac- 
tion, and antiferromagnetic superexchange, while ferro- 
magnetism and metallicity are favored by the double- 
exchange mechanismjii^ By simply changing the doping of 
a manganite you can go from ferromagnetic and metallic 
phases to antiferromagnetic and insulating. Phase sepa- 
ration is also commonly observed] "'^^d^d^ AU-manganite 
heterostructures are therefore a very interesting ground 
for analyzing the interplay between different phases, 
while minimizing the creation of defects at the interfaces 
thanks to the similarity in the lattice structure and chem- 
ical composition of the layers 

Motivated by the 'electronic reconstruction' effects ob- 
served in different oxide heterostructures and by the 
range of different phases that manganites show de- 
pending on their doping level, we have studied an 
all-manganite heterostructure formed by two different 
manganite insulators: Lao.sCao.sMnOa (LCMO) and 
CaMnOa (CMO). LCMO is a CE-type antiferromagnet 
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(consisting of ferromagnetic zig-zag chains in the x-y 
plane coupled antiferromagnetically to each other) and 
insulator in which the antiferromagnetic (AF) superex- 
change and the Jahn- Teller interaction open a gap in the 
band structure. On the other hand, CMO is a parent 
compound with no Cg electrons on the Mn sites and, 
therefore, a band insulator. CMO presents AF (tt, tt, tt) 
ordering (G-type) due to the superexchange between the 
t2g electrons. 

By looking at the band structure and the in-plane con- 
ductance of these heterostructures, we find that a ferro- 
magnetic (FM) and conducting interface can be formed 
for a range of reasonable parameters. We also find 
that by changing the electron concentration n in LCMO 
around half-doping, n — 0.5±(5, the charge in the metallic 
interface and the in-plane conductance can be manipu- 
lated. 

This paper is organized as follows. In Sec. |TT] we de- 
scribe the model we use to describe the manganite het- 
erostructures. The phase diagram, the characterization 
of the heterostructure ground state configuration, and 
the effect of doping around n = 0.5 are shown in Sec. lIIII 
Finally, we conclude in Sec. IIVI 
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II. MODEL 

We study the heterostructure depicted in Fig. [T] con- 
sisting of alternating CMO and LCMO layers a few unit 
cells wide. The background charge (namely, the charge 
on the AO planes) is in the CMO planes, +0.5 in the 
LCMO planes and, to achieve a symmetric distribution 
of the charge, we use the average -1-0.25 at the inter- 
faces. To study the phase diagram, different magnetic 
configurations on the CMO planes (either FM or G-type 
AF) are considered. On the other hand, the magnetic 
order in the LCMO planes is always considered to be of 
the CE-type: FM zig-zag chains in the x-y plane cou- 
pled antiferromagnetically with neighboring chains. As 
discussed below, FM planes in the LCMO layer are not 
favored energetically. 

The hamiltoniani^i^iiS^ includes the kinetic energy, 
a nearest neighbor antiferromagnetic coupling, and a 
Hartree term that takes into account the long range 
Coulomb interaction between all the charges in the sys- 
tem 

H — — ^ ^ fij^j,^'^i,-y^jn' ^ ^ "^AF^i^j -^Hartree j 

(1) 

where Cj^ creates an electron on the Mn i-site, in the 
eg orbital 7 (7 = 1, 2 with I = \x'^ - y^) and 2 = \3z^ - 
r^)). The hopping amplitude fij depends on the Mn 
core spins orientation given by the angles 9 and ijj via 
the double-exchange mechanism 

= cos(6»,/2)cos(6lj/2) + exp[j(V'>-Vj)]sin(6l,/2)sin(6l,/2), 

(2) 
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x,y, and z refer to the direction in the lattice. All the pa- 
rameters are given in units of t which is estimated to be 
~ 0.2 — 0.5 eV. Jaf is, in general, an effective antiferro- 
magnetic coupling between first neighbor Mn core spins 
which is different in the CMO and LCMO layers. For 
CMO, the antiferromagnetic coupling is a pure superex- 
change between the localized t2g spins, while for LCMO 
it would effectively include the interaction between the 
eg electrons and the latticci^ -ffHaitioc takes the form 
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(3) 

with Ri the position of the Mn ions, (n.^) — J^-yi^I j^i.-y) 
the occupation number on the Mn i-site, eZi the charge of 
the A-cation located at R^, and e the dielectric constant 
of the material. The relative strength of the Coulomb 
interaction is given by the parameter a = /aet [l8l |. 

We are assuming that the Cg Mn levels in the two dif- 
ferent layers are aligned, although a band offset due to 
a mismatch in the work functions of the two materials 
could occuri^ We expect our approximation to be rea- 
sonable as the two layers are chemically similar. 

We find the ground state configuration (at temperature 
T = 0) of the heterostructure by solving the hamiltonian 
in Eq. [T] self-consistently in a 4 x 4 x 12 system (6 unit 
cells of CMO and other 6 of LCMO, as depicted in Fig.[T]) 
with periodic boundary conditions in the three directions. 
Bulk behavior is always recovered at the center of the 
layers and the two CMO/LCMO interfaces (at planes 
i = 3 — 4 and planes i = 9 — 10) are independent of each 
other. 



III. RESULTS 

A. Phase diagram 

In Fig.[2]we show the phase diagram of the heterostruc- 
ture built by comparing the energies of four different con- 
figurations in the CMO layer. All the considered config- 
urations have a symmetry plane between planes i = 6 
and 7 in Fig. [T] (i) the configuration labeled OEM corre- 
sponds to G-type AF order on all CMO planes (j = 1 to 
3, and j = 10 to 12); (ii) IFM stands for one FM plane at 
each interface {i — 3 and 10); (iii) 2FM stands for two FM 
planes at each interface [i — 2 and 3 on one interface, and 
i = 10 and 11 on the other); and, finally, (iv) 3FM cor- 
responds to all FM CMO. We have checked that having 
FM planes in the CMO layer is more favorable energet- 
ically than having FM planes in the LCMO layer. This 
is mainly due to the fact that the Jahn- Teller coupling 
occurring in bulk LCMO, and absent in bulk CMO (with 
no Cg electrons), implies a larger effective AF superex- 
change in the former^ making AF order more robust in 
LCMO than in CMO. 
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FIG. 2: (Color online) Phase diagram for n = 0.5 in terms of 
the two parameters a and Jaf- Results are very similar for 
different values of n. The label in each region stands for the 
number of FM planes in CMO at the interface with LCMO. 
The 3FM region corresponds to a FM-CE heterostructure. 
In a wide range of reasonable values of Jaf (see text), the 
configuration with one FM plane between G and CE (IFM) 
is stable. 




plane index 



FIG. 3: (Color online) (a) Electron density at each plane for 
n = 0.5 and a = 2 for the configurations with no FM plane 
at the interface (OEM), with one FM plane at the interface 
(IFM), and with two FM planes (2FM). The curves for 3FM 
are very similar to 2FM. This illustrates how the electronic 
charge spreads across the interface. This spreading is limited 
by the Hartree potential shown in (b). 



As the LCMO layer has a fixed CE configuration in all 
planes, the phase diagram is independent of the J^™*-* 
coupling and it only depends on J^f"^ which is a pure 
AF superexchange {J^f" = Jaf 1 - 10 meV [ll). 
An estimation of Jaf in units of the hopping parameter 
t = 0.1- 0.5 eV would be Jaf ~ 0.02 - 0.1 1. In Fig. [H 
we plot the phase diagram as a function of a and Jaf- 
A likely value for a is in the range ~ 1 — 2 [l8l |. 

In Fig. Owe show the electronic density per plane and 
the Hartree potential for the OFM, IFM, and 2FM con- 



figurations. The results are shown for a = 2 (the curves 
are qualitatively similar for different values of a ~ 1 ~ 2). 
Fig. [3] illustrates the charge spreading across the inter- 
face and how it is limited by the Coulomb interaction 
between all the charges in the system, given by the aver- 
age on each plane of the Hartree potential, which at site 
i takes the form 

Wee« = ^E(^^^-^-^) . (4) 

The bulk values for the electronic density are recovered in 
the center of both CMO and LCMO layers. The curves 
for the 3FM (not shown) are almost coincident to the 
ones corresponding to 2FM. 

The phase diagram (Fig. [5]) shows that configurations 
with ferromagnetic planes are the ground state for the 
lowest values of the Jaf- For reasonable values of this 
parameter, the ground state is the IFM configuration. 
From now on, we focus on this particular case. 



B. Analysis of the metallic interface: density of 
states, conductance, and 'shadow' order. 

We analyze now the properties of the ground state IFM 
configuration by looking at the density of states (DOS), 
the in-plane conductance, and the charge distribution on 
each of the planes. As a reference, let us first discuss the 
density of states of a 2-dimensional CE system illustrated 
in Fig. m (a). For n — 0.5 the first band is completely 
full and the system is an insulator. For n > 0.5, the 
extra electrons would have to go to the non-dispersive 
band at energy = 0, but this requires a large energy, 
of the order of the gap, so the system would probably 
tend to phase separate with FM regions (this is the case 
of Pr„Cai_„Mn03 This phase separated system 

would also be insulating. For n < 0.5, assuming the CE 
order is not altered, the first band would be partially full, 
there would be a finite DOS at the Fermi energy, and the 
system would be conducting. However, the experimental 
observation is somewhat different: for n < 0.5 it is pos- 
sible that an incommensurate charge-orbital order arises 
opening a gap at the Fermi energy t^^'^^ leading to insu- 
lating behavior. The DOS of a 2-dimensional FM with 
two Cg orbitals has the peculiar shape shown in Fig. [3] 
(b). 

In Fig. [5] we show the DOS on planes i = 1 to 6 of 
the CMO/LCMO heterostructure in the IFM configura- 
tion. The Fermi energy (Ep) is signaled by the dashed 
vertical line crossing the six panels. The shape of the 
DOS of a perfect CE (as in Fig. 2] (a)) is recovered in 
planes i = 5 and 6, except for an overall shift due to the 
Hartree potential [Fig. [3] (b)]. The electronic density in 
plane i = 6 is exactly 0.5 (see Fig. [3]), the Ep is at the 
gap, and the plane is an insulator. On the other hand, the 
electronic density in plane i = 5 is slightly smaller than 
0.5, there is a finite DOS at Ep and, therefore, this CE 
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FIG. 4: (Color online) Density of states (DOS) for a perfect 
CE antiferromagnet (a) and a 2-dimensional FM (b). 
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FIG. 5: (Color online) Density of states (DOS) on each plane 
of the heterostructure (planes 1 — 6) in the IFM configuration 
for n = 0.5 and a = 2. The vertical lines signal the position 
of the Fermi energy. In the center of the LCMO layer (i — 
5, 6) we recover the bulk CE phase [see Fig |4ja)] which is 
an insulator, while the CE planes closer to the interface have 
states at the Fermi energy. The DOS is discrete, calculated 
with a large number of k-points, and the curves are smoothed 
out by giving a gaussian weight to the eigenvalues. 



plane conducts. The CE plane adjacent to the FM plane 
{i = 4) has a DOS very different from that of a perfect 
CE, in particular, the gap at the E^ closes and the plane 
conducts. The FM plane (i — 3) shows no gaps, as ex- 
pected. We have also performed numerical calculations 
via the Kubo formula^ii^ of the in-plane conductance 
that quantify these observations. In particular, by con- 
necting the planes one by one to the leads, we have cal- 
culated each plane's conductance. The results are shown 
in Fig. [S] (a) and show that three planes {i = 3, 4, and 
5) conduct. However, following the discussion in the pre- 
vious paragraph, it is clear that the finite conductance 
we find in the CE planes is an artifice of our approxi- 
mation of having fixed the CE magnetic ordering on the 
LCMO planes: it is possible (and consistent with experi- 
mental observations) that other complicate magnetic and 
charge-orbital orders or phase separation arose, leading 
to an insulating behavior. Exploring these other possibil- 
ities is out of our computing capabilities and we expect 
that imposing an insulating character on the doped CE 
phases would not affect the electrostatics of the system. 
The most reasonable interpretation of our results is illus- 
trated in Fig. [5] (b) where only the finite conductance in 
the FM plane is kept. 

Another important observation is that the DOS of the 
FM plane at the interface (i = 3) is very different from 
that of an isolated 2-dimensional ferromagnet [see Fig. [3] 
(b)]. This is a manifestation of the influence produced by 
the adjacent G and CE planes. This strong influence is 
also revealed by the appearance of a charge modulation 
('shadow' order) on the FM plane in the heterostructure, 
illustrated in Fig. [T] Here we show both the magnetic 
ordering (G, FM, and CE) and a qualitative representa- 
tion of the charge and orbital configurations for planes 
i = 2, 3, 4, and 5. Black dots correspond to sites with 
more charge than in gray dots. The elongated dots rep- 
resent orbital order. The charge modulation, which is 
absent in plane i — 5, arises in planes 3 and 4 due to the 
fact that some sites have neighbors with a parallel spin 
in the z-direction while others have antiparallel neigh- 
bors. Double exchange mechanism implies that hopping 
is not allowed in the case of antiparallel spin and, con- 
sequently, the charge is larger at those sites with more 
parallel neighbors. 



C. Selective doping on the LCMO layer 

We now explore the effect of doping the LCMO layer to 
produce a small increase or decrease of electronic density 
S (electron or hole doping, respectively). The considered 
heterostructure is still the same as in Fig. [T] but now the 
charge on the center slab is n — 0.5 ± 6 [(0.5 ± S)/2 at 
the interfaces]. The phase diagrams for the cases S ^ 
are very similar to the one for n = 0.5 shown in Fig. [2l 
However, due to the asymmetrj«2^i2^ on the electron or 
hole doping around n = 0.5, manifest on the shape of 
the DOS in Fig.[l](a), we expect, and indeed observe, a 
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FIG. 6: (a) Conductance versus the plane index for n = 0.5. 
The finite conductance on the CE layers (i = 4, 5, 8, and 9) 
is caused by our approximation of fixing the CE ordering on 
those planes, (b) In fact, we only expect the FM plane to 
conduct (see Sec. IIIIBl for discussion). 



qualitatively different behavior for n > 0.5 and n < 0.5 
in terms of the charge distribution. 

In Fig. [S] we show the electron density and Hartree po- 
tential for n = 0.46, n = 0.5 and n = 0.54. For n = 0.54, 
the electrons on excess of 0.5 cannot be accommodated 
in the center of the LCMO layer due to the gap in the 
DOS [Fig. g] so they go to the LCMO planes close to 
the interface and to the FM plane. On the contrary, the 
defect of electrons for the case n — 0.46 easily accommo- 
dates at the center of the LCMO layer and there is an 
overall decrease of the electron density at all planes. Con- 
sequently, the electron density in the FM plane increases 
with increasing n. 

The in-plane conductance G in the FM plane as a func- 
tion of n around half-doping is shown in Fig. [5] for a = 2. 
For this particular value of a, the conductance increases 
with both electron and hole doping around n = 0.5. How- 
ever, the shape of the curve is different for different values 
of a. This variability is due to the complicated structure 
(with many maxima and minima) of the DOS in the FM 
plane (see Fig.[S]). The structure on the DOS arises from 
the interplay between the DOS of a perfect CE [Fig. [4] 
(a)], that of 2-dimensional FM [Fig. g] (b)], and that of a 
G plane (single peak at energy^ 0), shifted relatively to 
each other by the Hartree potential. Changing n and/or 
a produce a change in both the Fermi energy and the 
Hartree potential so the peaks on the DOS move with 
respect to each other. As the conductance G is a mea- 
sure of the DOS at E^?, it is not surprising that small 
changes in n and a can produce relatively large changes 
in G. 
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FIG. 7: (Color online) Qualitative illustration of the magnetic 
ordering and charge distribution in different xy planes of the 
heterostructure for n = 0.5 in the IFM configuration. Black 
and gray dots represent points with more and less charge re- 
spectively. Elongated dots represent orbital order. The mid- 
dle CE planes (i = 6-7, not shown) behave as in bulk: they 
have homogeneous charge (0.5 per site) and orbital ordering 
that distinguishes the bridge from the corner sites in the zig- 
zag chains. They are also orbital polarized — because 
hopping is completely suppressed in the z-direction. Plane 
5 also has homogeneous charge although < 0.5. Plane 4 de- 
velops charge modulation with more charge on the zig-zag 
chains parallel to the FM plane {i — 3). Plane 3 does not 
have charge homogeneity, although it is FM, due to the gain 
in kinetic energy with the sites on the parallel zig-zag chains 
in plane 4. There is also a trace of orbital ordering in this FM 
plane, which mirrors the orbital ordering on the adjacent CE 
plane. Plane 2 gets some charge at the sites with spin parallel 
to the neighboring FM plane (only on the — orbitals). 



c 

03 
"O 

c 
o 

1_ 

o 

03 
03 



03 

O 
Q. 

03 
03 

CD 
X 




1 23456789 10 11 12 
plane index 

FIG. 8: (Color online) (a) Electron density per plane in the 
IFM configuration for a = 1.75 and three different values of 
n. (b) Corresponding Hartree potential. Due to the form of 
the density of states of a perfect CE plane [see Fig|4ja)] it is 
very expensive to add electrons to n = 0.5 so the charge goes 
to the planes with an electron density < 0.5. 
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FIG. 9; (Color online) In-plane conductance of the FM plane 
versus n around half-doping for a = 2. The overall shape of 
this curve depends on the value of a. 

IV. CONCLUSIONS 

The possibility of growing good quality oxide het- 
erostructures is opening a new playground to study both 
fundamental aspects of strongly correlated systems and 
the competition between their different phases, and po- 
tential applications of the new phases that arise at the 
interfaces. Manganites, with their very rich phase dia- 
grams, are excellent candidates to analyze all these is- 



sues. 

Here we have studied an all-manganite heterostructure 
consisting of a band insulator (CaMnOa) and a half- filled 
ordered manganite (Lao.sCao.sMnOs) in which the insu- 
lating behavior is driven by strong correlations in the 
system. We use a model which has been proved suc- 
cessful in reproducing the phase diagram of bulk man- 
ganitesi^ The Coulomb interaction between the charges 
in the system, introduced at a Hartree level, produces a 
spreading of the charge across the interface. This leads 
to the formation of a metallic and ferromagnetic plane at 
the interface, similarly to what has been observed in het- 
erostructures of Mott and band insulators,^ for a wide 
range of physical parameters. We have characterized 
this interface by its density of states and the in-plane 
conductance. The ferromagnetic metal at the interface 
shows charge and orbital modulation due to the inter- 
action with the ordered neighboring planes. Finally, we 
have also observed that small electron and hole doping in 
the Lao.sCao.sMnOs layer can produce large changes in 
the conductance of the ferromagnetic metallic interface. 
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Spain) and M.J.C. Ramon y Cajal program (MEC, 
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